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Inverse lift: a signature of the elasticity of complex fluids? 
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Abstract 

To understand the mechanics of a complex fluid such as a foam we propose a model experiment 
(a bidimensional flow around an obstacle) for which an external sollicitation is applied, and a local 
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response is measured, simultaneously. We observe that an asymmetric obstacle (cambered airfoil 
profile) experiences a downards lift, opposite to the lift usually known (in a different context) in 
aerodynamics. Correlations of velocity, deformations and pressure fields yield a clear explanation 
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q ■ of this inverse lift, involving the elasticity of the foam. We argue that such an inverse lift is likely 



common to complex fluids with elasticity. 

PACS numbers: 82.70.Rr, 83.80.Iz, 47.50. +d 
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A liquid foam exhibits "complex" be- polymer solutions ((J Q, 0, 0)), and is easily 

haviour under stress: it is elastic for small de- observable. In particular, a foam with only 

formation, plastic for large deformation, and one bubble layer (so-called "two-dimensional 

is easy to image, and image 

O ■ 

rich mechanical behavior is used in many of analysis yields information on all the geomet- 

the foams' applications, including ore sepa- rical properties of the foam. 
c$ ' ration by flotation in mines, drilling and ex- 



12, 



traction in oil industry, and cleaning in con- Weperform a Stokes experiment Q, 
fined media such as pipes A foam is a UJ, L3J, uA UM, i.e. we study the flow of 



convenient model to study constitutive rela- foam around obstac 
tions, since the microscale is the scale of bub- described in Ref. 



ues, using a set-up fully 
12j. Briefly, a tank is 



bles (not of molecules, as in most complex filled with a bulk solution obtained by adding 
fluids, such as emulsions jj, [J, colloids and 1% of commercial dish-washing liquid (Taci, 

Henkel) to desionised water. Its surface ten- 
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sion, measured with the oscillating bubble (zero-torque) stable equilibrium orientation 
method, is 7 = 26.1 ± 0.2 mN m _1 , and its by measuring the leading angle a, defined as 
kinematic viscosity, measured with a capil- the angle between the axis passing through 
lary viscosimeter, is 1.06 ±0.04 mm s -2 . Ni- the points x = ±1, y — in the Joukovski 
trogen is blown in the solution through a noz- equation (see caption of Fig. Q), and the 
zle or a tube at a computer controlled flow flow direction. The leading angle (which de- 
rate. This generates a foam, constituted by pends on the location of the contact point, 
a horizontal monolayer of bubbles of average and hence is not generic) is small and neg- 
thickness hn = 3.5 mm, confined between the ative: it decreases from —1° to —4° in the 
bulk solution and a g .a 8S top plate Q. The stud.ed ra„ g e of flow rate (Fig. B . 
foam is monodisperse (bubble area at chan- Fig. El reports the zero-torque orienta- 
nel entrance: A = 16.0 ± 0.5 mm 2 ) and its tion, and corresponding drag and lift mea- 
fluid fraction is around 10% (the evaluation surements, versus the flow rate. It evidences 
of this quantity in such a setup will be de- a non-zero drag at vanishing flow rate, which 
tailed in future work). It flows around an is the force required to trigger a steady mo- 
obstacle placed at the middle of the channel, tion of the foam with respect to the obsta- 
The obstacle is linked to a fixed base through c le, and appears more as a solid-like property, 
an elastic fiber; we thus measure the force On the other hand, the drag is an increasing 
exerted by the flowing foam on the obstacle affine function of the flow rate, as expected 



(precision < 0.1 mN) by tracking the obstacle [l2| [l3j (and its value is almost as low as 

displacement from its position at rest, using that of a non-cambered airfoil |l2j): this is a 

a CCD camera which images the foam flow conse que nce of the fluid-like properties of the 

from above. The flow rate is 50 ml min -1 , foam jl2j. 

and the average velocity 2.7 mm s" 1 , except pig. H also shows a striking feature: the 

in Fig. |21 lift is directed downwards. This is opposite 

Here, the obstacle is a cambered airfoil to the lift which appears (in an entirel y d iffer- 

(Fig. HJ). Like every obstacle in a flow, ent physical regime) in aerodynamics [17]. To 

the airfoil experiences a streamwise force, the our knowledge, this is the first time that such 

drag; but owing to its asymmetry, it also an inverse lift is experimentally evidenced, 

feels a torque, and a spanwise force: the lift. Does it originate from solid- or liquid-like 

The obstacle is free to rotate around the con- property? 

tact point with the fiber. We quantify its As a first hint, we note that the lift hardly 
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FIG. 1: Top view of the airfoil with a small 
part of the flowing foam. The superimposed 
white line indicates its boundary. It is a 

n 

Joukovski profile [17l| . which equation of shape 
writes x(t) = {1.56 + [(0.168 + cost) 2 + (0.25 + 
sint) 2 ]- 1 }(0.168 + cost) and y(t) = {1.56 - 
[(0.168 + cos t) 2 + (0.25 + sin t) 2 ]" 1 }(0.25 + sin t) , 
with lengths in centimeters, the angle t rang- 
ing from — 7r to 7T. Its length from leading to 
trailing edge is 5.2 cm. The arrows indicate 
the direction of drag and lift (flow from left 
to right), and the white dot marks the contact 
point between the airfoil and the fiber. The dot- 
ted lines show where the quantities displayed in 
Fig. 0] are evaluated. A movie is available at 
www-lsp/link/mousses-f ilms .htm. 

increases with the flow rate. To understand 
its physical origin, we now turn to the effect 
of the obstacle on the foam flow. Using bub- 
bles as passive tracers, we perform local mea- 
surements of the velocity, area and deforma- 
tion fields (which correlate respectively with 
viscous, surface tension and pressure contri- 
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FIG. 2: Forces exerted by the flowing foam on 
the airfoil versus flow rate: drag (□) and lift (•). 
Insert Leading angle (°) spontaneously selected 
by the airfoil versus flow rate (ml min -1 ); the 
standard deviation of its time fluctuations are 
plotted as error bars. 

butions to the stress ^,18, 3)- Their time 
averages are plotted around the airfoil (Fig. 
0J), and along two horizontal lines, 1 cm above 
and 1 cm below the airfoil (Fig. HJ). 

The velocity field shows that the convex 
regions of the airfoil constrict the flow. At 
the trailing edge's cusp, the velocity field is 
regular. It does not exhibit singularity, nor 
any qualitative difference with aerodynamics, 
where the velocity at a sharp trailing edge is 
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20j). 



continuous (Kutta condition 

The 3D compressibility of bubble gas is 
generally neglected (foams compression mod- 
ulus, of order of atmospheric pressure, is typ- 
ically three orders of magnitude larger than 
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FIG. 3: Experimental measurements of veloc- 
ity, area and deformation fields around the airfoil 
across the whole channel (width: 10 cm). Each 
measurement results from an average over a rep- 
resentative volume element (square box of side 
1.1 cm) and over time (750 successive images in 
steady regime). Arrows: velocity of bubble cen- 
ters of mass. Background gray levels: bubble 
area, with 12% relative variation between the 
most (dark grey) and less (white) compressed 
bubbles: the mean area corresponds to the con- 
tour line between the two gray levels at the left 
side of the figure. Ellipses: texture tensor, the 
major axis representing the direction and magni- 
tude of maximal bubble elongation (an isotropic 
region would be represented by a circle). 
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their shear modulus [l|). But here, thanks to 
the bulk solution in contact with the bottom 
of the foam, an increase of pressure increases 
the height of the bubbles (which equalise 
their pressure with the hydrostatic pressure 



of the bulk solution). Its effect is a de- 
crease of the visible bubble area. The present 
foam thus has an effective 2D compressibil- 



ity, equal to (pgW 



2,9x KT 2 Pa 



-i 



[12]. 



For such a compressibility, bubble area vari- 
ations act as a passive tracer of the pressure 
field: they are large enough to be measur- 
able, and small enough not to perturb the 
flow. From area measurements, we can deter- 
mine the net contribution of pressure to the 
force: it writes |l 21] Fv = -pgA\h\ § din/ 'A 2 , 
where p = 1.00 x 10 3 kg m -3 is the volumetric 
mass of the solution, g = 9.8 m s~ 2 the grav- 
ity acceleration, ho and A the average values 
of the bubbles' depth and area, the integral 
being taken over the contour of the airfoil (n 
is the outwards normal of the contour and 
d£ its length element). We have measured 
Fp. It contributes for 0.20 mN to the drag, 
and for 0.74 mN to the downwards lift. It is 
worth noting that the bubbles' pressure in- 
creases where the flow accelerates (Figs. El 
and |H), contrary to Newtonian fluids in iner- 

n 

tial flow 20 : this is a clear signature of foam 
elasticity. 

We quantify the deformations of bub- 
bles, visible on an image ^|, by measur- 
ing the_ (always symmetric) texture tensor 
defined as: M = (£ <g> 0) = 



((£ 2 X ), (£J y ), (£J y ), (£ 2 y )). It only requires to 
measure the bubble edge vectors £ linking two 
neighbouring vertices; the average is taken 
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FIG. 4: Experimental measurements of veloc- 
ity, area and deformation fields (a) 1 cm above 
and (b) 1 cm below the airfoil. Same data as 
Fig. |21 The streamwise component v x of the 
velocity (□), the bubble area (•) and the yy 
component of the texture tensor, M yy (A) are 
adimensioned by their value at the channel en- 
trance (2.7 mm-s" 1 , 16 mm 2 and 3.5 mm 2 , re- 
spectively), and we represent their variations rel- 
ative to these values versus the streamwise co- 
ordinate x relative to the leading edge. Vertical 
dots indicate the leading and trailing edges, 
inversion below the airfoil occurs 1.8 cm after the 



over a representative volume element. Fig. El 
shows the elongation of the bubbles on both 
sides of the airfoil: vertical stretching in the 
concave region of the airfoil (below the trail- 
ing edge); horizontal stretching in the con- 
vex regions (above the airfoil and below the 
leading edge). This qualitatively different 
behaviour appears clearly on the M yy plots 
(Fig. ij). 

Correlations visible on Fig. El and 0] yield 
a physical explanation of the downwards lift. 
In convex regions, above the airfoil or below 
aX x < 1.8 cm, the flow is constricted and 
accelerates, and bubbles stretch streamwise. 
Since the elastic stress in foams is due to sur- 

ated to the orientation 

3, EH 



face tension and corre 



of bubble edges 
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of main elastic stress is streamwise. In the 
concave region, below the airfoil at x > 1.8 
cm, the picture is reversed, and the direc- 
tion of main elastic stress is spanwise; this 
contributes to a downwards lift, like the bub- 
bles' pressure contribution described above. 
The net balance is a resulting downwards lift. 
This explanation is in principle only valid in 
a quasistatic regime, but since the lift does 
not increase significantly with the flow rate 
(Fig. 2), it seems to constitute the essen- 
tial ingredient in the studied range of flow 
rate. Moreover, the "Tl" bubble neighbour- 
swappings [1] saturate the maximum value of 
deformation but do not fundamentally affect 



this mechanism. The lift thus appears mainly sible other contributions. First, to exclude 

as an elastic effect, typical of a solid-like be- possible artifacts linked with the present set- 

haviour. up, namely bubble 3D geometry and their 

To estimate the elastic contribution to the effective 2D compressibility, S. J. Cox (pri- 

lift, we approximate the foam by a 2D one, vate communication) performed simulations 

with a line tension A = 2jh = 0.18 mN, prob- of a true 2D, incompressible foam flow using 



He unam- 



ably a strong underestimation (due to 3D ef- the Surface Evolver software 

fects, the effective height and thus line ten- biguously observed the same bubble stretch- 

sion is probably larger). The 2D elastic stress ing and downwards lift, due both to the bub- 

in such a 2D foam would write as: ble edges' surface tension and to the pres- 

<7 e i = Xp{£ where / denotes a bubble sure contribution. Second, the confinement 

edge, and p = 3/ A the areal density of edges of the foam by the sides of the channel is ex- 

(there are in average six edges per bubble pected to play a role: this is always the case 

each one being shared by two bubbles). We in 2D flows around obstacles^dther Newto- 

n nn 

24j or non-Newtonian |15j, [16|. How- 



estimate an elastic contribution of 0.05 mN man 

to the drag, and 0.23 mN to the lift. Adding ever, the relevant parameter is the logarithm 

pressure and elastic contributions indicates of the channel width to obstacle size ratio, 

values of 0.25 mN for the drag and 0.97 mN and experimental studies of the drag exerted 

for the downwards lift, of the same order but by a flowing foam on obstacles show weak 

lower that the measured forces. variations of the drag with the ratio obstacle 

At this stage, we can propose an expla- size/channel width |n|,|lj. Hence, we expect 
nation for the low dependence of the lift on a weak quantitative (and no qualitative) ef- 
the flow rate. Since the foam slips along the fect of the channel width on the lift. Third, 
airfoil, in the lubrication films between the the aerodynamic lift scales like U sin(a + (3) 
airfoil and the surrounding bubbles there ap- [li M, where U denotes the relative veloc- 
pears strong velocity gradients Q; they are ity of the flow and the obstacle, and p the 
perpendicular to the airfoil boundary, hence P urel y geometric camber angle. For our air- 
mainly perpendicular to the flow. The resul- foil this an S le ec i uals P = 14 - 5 °! hence > even 
tant of the viscous friction thus contributes if " is negative, a + (3 remains positive, and 
much more to the drag than to the lift. this cannot explain our observations. Fourth, 

To discuss whether this lift is a true ef- there is an avera § e Pressure gradient VP due 

fect of foam physics, we must examine pos- to the dissipation of flowing foam It 



equals 40 Pa m _1 for a flow rate of 50 ml 
min -1 , and reaches 1.7 x 10 2 Pa m -1 at the 
highest studied flow rate (565 ml min -1 , see 
Fig. |2J) . It slightly prestrains the bubbles be- 
fore arriving on the obstacle, but this does 
not qualitatively affect the main features of 
the deformation. It also adds an Archimedes 
thrust-like downstream contribution to the 
drag: II = Sh VP (S = 7.74 mm 2 : sur- 
face of the airfoil), which is not negligible 
(VP = 0.11 mN at 50 ml min -1 and up to 
0.46 mN at 565 ml min -1 ), but which does 
not contribute to the lift. 

How generic is the inverse lift? First, it 
is compatible with other phenomena (for in- 
stance, die swell, Weissenberg rod-climbing 
effect p Q, 0, a], sedimentation ofparti- 
cles |26j, or inverse Magnus effect [271 ]) ob- 
served or predicted with non-Newtonian flu- 
ids which act in the opposite sense to New- 
tonian fluids in inertial flow j^. More pre- 
cisely, wherever the pressure of a Newtonian 
fluid would push an obstacle, the normal 
stress in a viscoelastic fluid pulls it; for in- 
stance it can change from compression to ten- 
sion at the trailing edge j^, in agreement 
with what we observe here. Second, prelim- 
inary studies of the flow of a second-order 
fluid on the same airfoil profile show unam- 
biguously the inverse lift Note however 
that (contrary to the present case) for fluids 
with zero yield stress, the lift would be ex- 



pected to vanish at vanishing flow rate, and it 
may increase significantly with the flow rate 
if normal stress differences do (2^ . Third, we 
experimentally let an asymmetric object (a 
truncated portion of a disk, with a circular 
side and a straight one) settle under grav- 
ity in a model viscoelastic fluid (0.5% w:w 
cellulose solution) confined between vertical 
plates of glass. The object does feel a lift di- 
rected from the most to the less convex side. 
Several arguments thus suggest that such an 
inverse lift is expected to be generic to other 
fluids which can store elasticity. 
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